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6. Matrix integrals and counting planar diagrams 

6.1. The number of planeir gluings. Let us return to the setting of §4. Thus, we have a potential 
U{x) = — J2j>o9j^'' /j (with gj being formal parameters), and consider the matrix integral 

Jt)N 

Let ZN^h) = Zi\f{h/N). We have seen that 

where Woo is given by summation over planar fat graphs: 



In particular, the coefficient of Y\ g^* is (up to a power of H) the number of (orientation preserving) 
gluings of a fat graph of genus zero out of a collection of fat flowers containing rij i-valent flowers for 



I 'n 



each i, divided by Y[ 

On the other hand, one can compute Woo explicitly as a function of gi by reducing the matrix integral 
to an integral over eigenvalues, and then using a fundamental fact from the theory of random matrices: 
the existence of an asymptotic distribution of eigenvalues as TV ^ oo. This approach allows one to 
obtain simple closed formulas for the numbers of planar gluings, which are quite nontrivial and for 
which direct combinatorial proofs were discovered only very recently. 

To illustrate this method, we will restrict ourselves to the case of the potential U{x) = x'^/2 + gx^ 
(so gi = —Ag and other g^ = 0), and set fi = 1. Then Woo = X]n>i Cn{— g)"' / n\ , where c„ is a number 
of connected planar gluings of a set of n 4-valent flowers. In other words, c„ is the number of ways (up 
to isotopy) to connect n "crosses" in the 2-sphere so that all crosses are connected with each other, and 
the connecting lines do not intersect. 

Exercise. Check by drawing pictures that c\ = 2, C2 = 36. 

Theorem 6.1. (Brezin, Itzykson, Parisi, Zuber, 1978). One has 

c„ = (12)"(2n-l)!/(n + 2)! 

6.2. Proof of Theorem 6.1. Let us present the proof of this theorem (with some omissions), We will 
assume that g is & positive real number, and compute the function Woo (5) explicitly. We follow the 
paper of Brezin, Itzykson, Parisi, and Zuber "Planar diagrams". Comm. Math. Phys. 59, p. 35-51, 
1978. 

The relevant matrix integral has the form 



Passing to eigenvalues, we get 



^ JN{g) 



where 



Jiv(O) ' 

Jn{9)= I e-^(^^?/2+«^^^)n(A.-A,)^dA. 



Thus, Woo(ff) = E{g) - E{Q), where E{g) = lim^v^oo A^"^ lnJjv(5)- 

Proposition 6.2. (steepest descent principle) E{g) equals the maximal value of the logarithm of the 
integrand. 

The proposition says, essentially, that the integrand has a sufficiently sharp maximum, so that the 
leading behavior of the integral can be computed by the steepest descent formula. We note that we 
cannot apply the steepest descent formula without explanations, since the integral is over a space whose 
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dimension grows as the perturbation parameter 1/A'' goes to 0. In other words, it is necessary to do 
some estimates which we will omit. We will just mention that for .9 = 0, this result can be derived from 
the explicit evaluation of the integral using Hermite polynomials (see §4). For the general case, we refer 
the reader to the book of P. Deift "Orthogonal polynomials and random matrices: a Riemann-Hilbert 
approach" . 

The integrand K{Xi, . . . , Aat) = e^^^^ ]^^^^.(Aj — A^)^ has a unique maximum, because 

it is logarithmically concave (check it!). The maximum of the integrand is found by equating the partial 
derivatives to zero. This yields 

(16) El^4^ = ^(^^^ + 25A?). 

Let Ai < A2 < • • • < Ajv be the unique (up to permutations) solution of this system of equations. 

Proposition 6.3. The normalized counting measures -j^ ^ 5{x — Xi) converge weakly to a measure 
IJl{x) = f{x,g)dx, where f{x,g) is a continuous function, supported on a finite interval [—2a, 2a], and 
differentiable on this interval. 

For the proof we again refer the reader to the book of P. Deift (p. 132 and later). We note that for 

g = 0, hy Wigner's semicircular law, a = 1 and f{x,0) = ^ 1/4 — ; so f{x,g) = ^\/4 — + 0{g). 

Now our job will be to find the function f{x,g). Passing to the limit in equation 16 (which requires 
justification that we will omit), we get 

/2a 1 
/(x, g)dx =-y + 2gy\ \y\ < 2a 
-2a y ~ ^ ^ 

where the integral is understood in the sense of principal value. 

This is a linear integral equation on f(x,g), which can be solved in a standard way. Namely, one 
considers the analytic function F{y) = /_2a ':^zr^f{x,g)dx for y in the complex plane but outside of the 
interval [—2a, 2a]. For y G [—2a, 2a], let F+{y), F-{y) denote the limits of F{y) from above and below. 
Then by the Plemelj formula, the integral equation implies 

^-{F+{y) + F_{y))=^-y + 2gy\ 

On the other hand, F+{y) = F-{y). Hence, YieF+iy) — Re(i^_(y)) = ^y + 2gy^. 

Now set y = a{z + z~^). Then, as y runs through the exterior of [—2a, 2a], z runs through the 
exterior of the unit circle. So the function G{z) := F{y) is analytic on the outside of the unit circle, 
with decay at infinity, and I{jeG{z) = \a{z + z^^) + 2ga^{z + z~^)^, \z\ = 1. This implies that G{z) is 
twice the sum of all negative degree terms of this Laurent polynomial. In other words, we have 

G{z) = Aga^z-^ + (a + I2ga^)z-'^. 

This yields 

Now f{y, g) is found as the jump of F: 

f{y,9) = -il+ 45«' + 2<72/') V4a2 - y^. 

It remains to find a in terms of g. We have yF{y) ^ 1, y — > 00 (as / f(x,g)dx = 1), hence 
zG{z) 1/a, z ^ 00. This yields 1/a = a + I2ga'^ , or 

125a'* + a^ - 1 = 0. 

This allows one to determine a uniquely: 

/ (l + 48g)V^-l 
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Now let US calculate E{g) . It follows from the above that 

/2a p2a /'2a -j^ 

/ lT^\x-y\fix,g)f{y,g)dxdy- {-x^ + gx^)f{x,g)dx. 
-2a J -2a J -2a ^ 

On the other hand, let us integrate the integral equation defining /(x, g) with respect to y (from to 
u). Then we get 

2a 

(In |a; — wj — In |x|)/(x, g)dx = —u^ + gu^. 

-2a 2 

Substituting this into the expression for E{g), wc get 

{\n\u\- -u^ - -gu^)f{u,g)du 

-2a 4 ^ 

Since f{u, g) is known, this integral can be computed. In fact, can be expressed via elementary functions, 
and after calculations we get 

i;(5)-£;(0) = lna--L(a2-l)(9-a2). 
Substituting here the expression for a, after a calculation one finally gets: 

^(.,)-E(0, = |(-12,)^i|-|. 

This implies the required formula for c„. 



